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Abstract 

We discuss two lagrangian interacting dark energy models in the context of the holographic 
principle. The potentials of the interacting fields are constructed. The models are compared 
with CMB distance information, baryonic acoustic oscilations, lookback time and the Constitution 
supernovae sample. For both models the results are consistent with a non vanishing interaction in 
the dark sector of the universe - with more than three standard deviations of confidence for one of 
them. Moreover, in both cases, the sign of coupling is consistent with dark energy decaying into 
dark matter, alleviating the coincidence problem. 
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I. INTRODUCTION 



In the last years, there have been several papers where an interaction in the dark sector of 
the universe is considered . A motivation to considering the interaction is that dark 

energy and dark matter will evolve coupled to each other, alleviating the coincidence problem 
A further motivation is that, assuming dark energy to be a field, it would be more natural 
that it couples with the remaining fields of the theory, in particular with dark matter, as it 
is quite a general fact that different fields generally couple. In other words, it is reasonable 
to assume that there is no symmetry preventing such a coupling between dark energy and 
dark matter fields. Using a combination of several observational datasets, as supernovae 
data, CMB shift parameter, BAO, etc., it has been found that the coupling constant is 
small but non vanishing within at least la confidence level [l] - J4]. In two recent works, 
the effect of an interaction between dark energy and dark matter on the dynamics of galaxy 
clusters was investigated through the Layser-Irvine equation, the relativistic equivalent of 
virial theorem Q. Using galaxy cluster data, it has been shown that a non vanishing 
interaction is preferred to describe the data within several standard deviations js]. However, 
in most of these papers, the interaction term in the equation of motion is derived from 
phenomenological arguments. It is interesting to obtain the interaction term from a field 
theory. Some works have already taken a step in such a direction {4], jg], (9]. On the other 
hand, scalar fields have been largely used as candidates to dark energy. They naturally 
arise in particle physics and string theory. Good reviews on this subject can be found in 



10| . A motivation to use scalar fields as candidates to dark energy is that their pressure 



can be negative, making possible to reproduce the recent period of accelerated expansion of 
the universe. For example, for canonical scalar field, the equation of state parameter varies 
between 1 and —1. Dark energy modeled as a canonical scalar field is called quintessence 
and was investigated, for example, in jj, For tachyon scalar field, the equation of state 
is always negative. The tachyon field has been studied in recent years in the context of string 
theory, as a low energy effective theory of D-branes and open strings 12[. Tachyon field as 
dark energy was studied, for example, in J], [3], [ijj. The first question about scalar fields 
concerns the choice of the potential. Common choices are the power law and the exponential 
potentials. However, these choices are in fact arbitrary. In principle, any other form for the 
potential which leads to recent accelerated expansion would be acceptable. 



On the other hand, it is possible that a complete understanding of the nature of dark 
energy will only be possible within a quantum gravity theory context. Although results for 
quantum gravity are still missing, or at least premature, it is possible to introduce, phe- 
nomenologically, some of its principles in a model of dark energy. Recently combinations of 
quintessence, quintom and tachyon models with holographic dark energy had been proposed 
-in 15 1, 16] and [yj], respectively. Specifically, by imposing that the energy density of 
the scalar field must match the holographic dark energy density, namely pa = Sc^M^L^ 2 , 
where c is a numerical constant and L is the infrared cutoff, it was demonstrated that the 
equation of motion of fields for the non-interacting case reproduces the equation of motion 
for holographic dark energy. In fact, to impose that the energy density of scalar field must 
match the holographic dark energy density corresponds to specify its potential. This can 
be seen as a physical criterion to choose the potential. Here, we generalize this idea for two 
kinds of interacting scalar fields. 

II. THE MODELS 

We consider the general action 

s = j d 4 xV=g l^R + £ v (x) + -[^V^ - f^y*] - (m - /fy>)M + ]T L 3 (x) 

(i) 

where Mpi = (8nG) 1 is the reduced Planck mass, R is the curvature scalar, C v (x) is, 
unless of the coupling term, the lagrangian density for the scalar field, which we will identify 
with dark energy, \I> is a massive fermionic field, which we will identify with dark matter, 

(3 is the dimensionless coupling constant and Cj (x) contains the lagrangian densities for 

j 

the ramaining fields. Notice that, in this work, we will only consider an interaction of dark 
energy with dark matter. If there was a coupling between the scalar field and baryonic 
matter, the corresponding coupling constant (3b should satisfy the solar system constraint 



[17| 



0b < 10- 2 . (2) 

We assume = 0, which trivially satisfy the constraint (T5]). 

We consider two kinds of scalar fields: the canonical scalar field, or quintessence field, for 
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which 

L 9 (x) = l -d^y ~ V(<p) , (3) 
and the tachyon scalar field, for which 

C v (x) = -V{<p) y/1 - ad^d^ , (4) 

where a is a constant with dimension MeV~ 4 . Notice that in both cases, we assume a 
Yukawa coupling with the dark matter field 



A. Quintessence field 

For the quintessence field, C v (x) in the action (pQ) is given by ([3]). From a variational 
principle, we obtain 



iY*V M tf - M*^ = , (5) 

+ M*% = , (6) 

where M* = M - f3ip, and 

V^ + ^M^W. (7) 

Eqs. (jSJ) and (jSJ) are, respectively, the covariant Dirac equation and its adjoint, 
in the case of a non vanishing interaction between the Dirac field and the scalar 
field (p. For homogeneous fields and adopting the Friedmann-Robertson- Walker metric, 
(7^=diag(l, —a 2 (t) , —a 2 (t) , —a 2 (£)), where a 2 (t) is the scale factor, eqs. jSJ) and © lead 
to 

d{a 3 ^) 



dt 



which is equivalent to 



and (JTj) reduces to 



f^ = f ^o(^) 3 (8) 



(p + 3Hp + = /3M , (9) 



where if = - is the Hubble parameter. 
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From the energy-momentum tensor, we get 

= \tf + V(<p) , (10) 
P<p = 1<P-V(<p), (11) 
p* = M*M , (12) 

= . 

From (fTUj) and (TTTj) we have w v = ^ = |^ 2+ ^|^ . Deriving ( ITOj) and (|T2|) with respect to 
time and using (JBJ) and ©, we obtain 

p v + 3#p^ + 1) = /3<pt # (-) 3 (13) 

and 

+ 3#p* = -/3^Wo (^) 3 , (14) 

where the dot represents derivative with respect to time. 
For baryonic matter and radiation, we have respectively 

p b + 3Hp b = (15) 

and 

p r + 3Hp r (u r + 1) = , (16) 

where u r — |. Eqs. f lT5|) and ffTB"]) implies that p b = and p r = respectively. The 
subscript denotes the quantities today. We are considering the radiation as composed 
by photons and massless neutrinos, so that p r o = (1 + 0.2271iV e //) p 70 , where N e ff = 3.04 
is the effective number of relativistic degrees of freedom and p 7 o is the energy density of 
photons, given by p 7 o = j^Tcmb, being Tqmb = 2.725-K" the CMB temperature today. The 
Friedmann equation for a flat universe reads 

a / 2 a 3 a 4 



if 2 



3M 2 , 



V a / 2 cr a 4 



(17) 



In order to determine the dynamics of the interacting quintessence field, it is necessary 
to specify the potential V((p). Instead of choosing an explicit form for V(ip), we will specify 
it implicitly, by imposing that the energy density of the quintessence field, given by ffTO]) . 
must match the holographic dark energy density, p\ = 3c 2 Mj^L -2 , where c is a numerical 
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constant and L is the infrared cutoff. The evolution of the interacting quintessence field with 
redshift will be given by the equation of evolution for the holographic dark energy density, 
with a certain expression for the equation of state parameter u v . In fact, we will see that 
imposing the energy density of the quintessence field to match the holographic dark energy 
densit y le ads to an expression for the potential. 

it has been argued that, in order that holographic dark energy drives the recent 



In 



1* 



period of accelerated expansion, the IR cutoff L must be the event horizon Rh = a (t) J t °° ^pj. 
Substituting Rh in the expression of the holographic dark energy, we get Rh = — 



therefore, 



dt' 



a (?) a (t) ' 

Differentiating both sides with respect to time, using the Friedmann equation (1171) together 
with conservation equations ( TT3T) and (fill) , we obtain 



dz 



a. 



l + z 



+ svl^u^ + n r + 1 



(18) 



Equation (1T81) is just the equation of evolution for the holographic dark energy (18J. Using 
the Friedmann equation (fTTj) . the conservation equations (fT5l) and (|T6l) can be written as 



dz 



l + z 



(3^^ + n r 



and 



dfl r 
dz 



[3Q v u v + Q r — 1) 



(19) 



(20) 



We define r = y-. Deriving r with respect to time, using ( 113"]) . flHJ) , ( JTUj) and (fTTj) we 



obtain 



r = SHrojui — sign 



0(l + r) /l + uv- 



VSM Pl H V fi v 



l + z 
l + z 



(21) 



We can rewrite ^o^o in terms of observable quantities. In fact, by imposing that the dark 



matter density today matches the observed value, we obtain ^/q^o 



M-p<f> 

The sign of (ft is arbitrary, as it can be modified by redefinitions of the field, <p — > —ip, and 
of the coupling constant, (5 —> —{3. Noticing that r = i - ^-^-^ ^ wg can SUDS titute r anc [ 



r in (1211) by O^, f2 r , fi^, Q b and fl r . Using ( ITHj) . ( JT9j) and (120]) we obtain, after some 
algebra 



00<p (z) 



1 

3 



"3c 



\ 



7 (z) 2 + 4 1- 



(22) 
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where 



with 



1 S7i/r (1 — ^<£>0 — ^60 — ^ro) fl + Z 



j(z) = -^5M P r yU ""^ , (23) 



E( )- H ^ / [(I ~ SAcp) (1 - ggp - n w - n r0 ) + n b0 ] n +z \ 3 Q r0 ( 1 + z V 

[z) ~ H V \i + «J i-^Vi + W 

(24) 

where Ay? (2;) = tp (z) — (po and 5 = M _f^ yo is an effective coupling constant. Notice that, if 
5 = 0, (|22|) reproduces the equation of state parameter obtained in fis| . 
The evolution of the quintessence field is given by 

dp y/3Mpiy/n ip (z)(l + u lf> (z)) 

Tz = TTz (25) 

From ffT8]) . ffT9~]) . f[20|) and fl25|) we can calculate the evolution with redshift of all observ- 
ables in the model. If we wish to calculate the time dependence, we need to integrate the 
Friedmann equation ([17]) . which can be written in the form 

dt 1 



dz H E(z)(l + z) 
From f lTOl) . we can compute the potential V (z) as 

PcO * 

where p c0 = 3M pi Hq, E (z) is given by (1241) . oj v (z) is given by fl22l and (z) is the solution 
of ( fl8l) . From ( 1261) and ( 1251) . we can compute ^(<£>). 

Here it is worth saying that in the holographic dark energy model, in the non interacting 
case - (1221) with 5 = - can be less than —1. However, as already mentioned in [ijj], if 
we wish that the holographic dark energy is the quintessence field, then because (1251) . oj v 
must be more than —1. Nevertheless, in the interacting case considered here, due to the 
fact that Up depends explicitly on ip, u v can not be less than —1. On the other hand, the 
square root in ( 1221) must be real. We can verify that is real and > — 1 if (i) < 1 



or (ii) > 1 and \S\ Mpi > ^n^n^-n - \j — 1. However, case (ii) is irrelevant, 

as it corresponds to large values oi\8\ Mp\. For example, if Q v o = 0.7 and c = 0.8, we have 
\5\ Mpi > 1.69. Below, we will see that the observational data constrain \5\ Mpi ~ 10 _1 . In 
order that be real for all future times, as Cl v —> 1, it is necessary that c > 1. 
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In(a) ln(a) 

FIG. 1: (a) Equation of state parameter of the holographic quintessence model, for c = 0.85 and 
5Mpi = —0.1 (red dashed line), 5Mpi = (black solid line) and SMpi = +0.1 (green dotted line), 
(b) Full range of the equation of state parameter for 5Mpi = +0.1. 

It is interesting to notice that the condition y ^^ p0 < 1 is precisely the same one for which 
the entropy of the universe increases [18|]. As f2 v — > 1 in the future, it is necessary that 
c > 1. Therefore, the condition for u v be real is precisely the same one for the entropy to 
increase. So the model respects the second law of thermodynamics. 

In figure [1] we see the evolution of the equation of state parameter u v with the scale factor 
a. For the non interacting case, 5 = 0, we have u v — > —1/3, as Q v (z) <C 1 for z 1. For 
5 < 0, u^ > — 1 in the matter era, then approaches —1 in the radiation era. For 5 > 0, u^, will 
eventually turns out positive and possibly u v 3> 1 in very early times, as in the case showed 
in figure [U This behaviour is explained as follows. In the matter era, E 2 (z) ~ (1 + z) 3 so 
that j(z) ~ I — From ( !T8|) < 0, so \j(z) \ increases with redshift. This increasing of 
|7 (z)\ will continue until the radiation era, when E 2 (z) ~ (1 + zf and Vt v (z) ~ (1 + z)~ 2 
so that |7(^)| — > cte. Tipically this constant will be much more than one. Therefore, for 



high redshifts (z) ~ — | + | 7 + bl + ^ 



. If S < then 7 < and tu^ (z) ~ —1 in the 



radiation era. If 5 > then 7 > and co^ (z) ~ | + |7 2 — > cte. Notice from fT2o]) that if 
Up (z) > 1 then V < 0. In order to avoid it, we impose the condition u v < 1 for all z. This 
condition is satisfied if 5M P i < jz^~z^~\/^ip- As we have Vl^ <ti 1 as z increase, there will 
be an abrupt decrease in the positive tail of the probability distribuction of 5, as well as in 
the confidence regions of S with other parameters. 

In figure [2], V((p) is shown for some values of 6 and c. Notice that there is a region where 
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FIG. 2: Potential of the holographic quintessence field V(<p), in units of p c o = 3M pi Hq. ip is in 
units of Mpi. The solid lines are for c = 0.95, the dashed ones are for c = 1.1 and the dotted are 
for c = 1.25. For each value of c the curves from right to left are for 5Mpi = — 0.1 (red), SMpi = 
(black) and 5Mpi = +0.1 (green), respectively. 

V(p) is almost constant, that is, there is a slow-roll region. As we chose ip positive, then 
ip evolves to this slow-roll region. However, if we had chosen p> negative, then because the 
right hand side of ( )25|) would have the opposite sign, so dV ^ would have also the opposite 
sign and again <p would evolve to the slow-roll region. Notice also that for SMpi = +0.1, 
the potential is negative in the past. 

The equation for evolution of p> (1251) can be written in an integral form as 



Since the model depends on Aip - through E (z) - and neither on ip nor on po, then it is 
independent of <pq. In other words, (po is not a parameter of the model and can be chosen 
arbitrarily. Therefore, the parameters of the model are S, c, h, f2 60 and Q^q. 

B. Tachyon scalar field 

In the case of dark energy modeled as the tachyon scalar field, (x) in the action (CQ) is 
given by (jl]). From a variational principle, we obtain 




i^V^ - M*^ = 



(27) 



+ m*^ = o , 



(28) 
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where M* = M — /3ip, and 

&<p{VM&><p l dlnV(ip) flM r _ — 

V ^ + a + a—d^ = ^)^* • (29) 

The equations of motion for \l/ and ^ (J27j) and (|28|) are the interacting covariant Dirac 
equation and its adjoint, respectively, i. e., ff27|) and f )28|) are almost the same as eqs. 
([5]) and ([6]), the only difference is that the scalar field <p in M* now is the tachyon 
field. For homogeneous fields and adopting the Friedmann- Robert son- Walker metric, 
(?^=diag(l, —a 2 (t) , —a 2 (t) , —a 2 (t)), where a 2 (t) is the scale factor, ( )29l) reduces to 



= -(! - "V? 2 ) 



1 dlnV(<f) . 



3//<p yl — a<p 



(30) 



a d(f ' aV(ip) 

whereas for the fermions, the equations of motion will reduce to eq. ([8]) , as already obtained 



above: 



^=^ (^y . (8) 



From the energy-momentum tensor, we get 

P v = f. = , (31) 
y 1 — a(f z 

P^ = -V{ V Wl-a<p 2 , (32) 
p* = M*M , 
Pa, = . 

From ( 13~TT) and ( 1321) we have oj v = ^ = a<p 2 — 1. Deriving ( 13~T]) and (1521 with respect to 
time and using (1301 and (JE]), we get 



Ptp + 3Fp^(uv + 1) = /9£* *o (^) 3 (33) 



and 

p* + 3#p* = -/9ptf *o (^) 3 , (34) 

where the dot represents derivative with respect to time. 

For baryonic matter and radiation, the conservation equations are the same as in the 
quintessence model. We have 

Pb + 3Hp b = (15) 
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and 



p r + 3Hp r (oj r + 1) = 



PrO 



(16) 

respectively. The 



where u r — |. Eqs. f lT5|) and f|T6|) implies that pb = and p r 
subscript denotes the quantities today. We are considering the radiation as composed 
by photons and massless neutrinos, so that p r o = (1 + 0.2271iV e //) p 70 , where N e ff = 3.04 
is the effective number of relativistic degrees of freedom and p 7 o is the energy density of 

2 

photons, given by p 7 o = j^Tcmb, being Tcmb = 2.725K the CMB temperature today. The 
Friedmann equation for a flat universe reads 



H 2 



3M 2 pi 



+ 



P60 PrO 

1 Q /I 



(35) 



a 1 a/I - oap 2 

Now, as done in the case of the quintessence field, we will identify the tachyon energy 
density ( 13~TT) with the holographic dark energy density p\ = 3c 2 M Pl L~ 2 . From a similar 
reasoning, we obtain again the equations f[T8"j) . f lT9"j) and (12"0~|). 







[2^ 




1 + z 






drib 


Sib 




dz 


1 + z 



+ 31^^ + Vl r + 1 



and 



Also we obtain 



3Hru,, 



dQ r 

dz 



sign 



l + z 



(3fi^ + Q r ) 



(3^c^ + £l r — 1) . 



/3(l + r)^yr+^ = 



3M Pl ^H 2 



l + z 
l + ^o 



(18) 
(19) 

(20) 
(36) 



with r = y-. The sign of <p is arbitrary, as it can be modified by redefinitions of the field, 
if — > — (f, and of the coupling constant, (3 — > —(3. We can rewrite ^0^0 m terms of observable 
quantities, by imposing that the dark matter density today matches the observed value. We 
obtain M^q^o — 3Mp '" H °( 1 ^ b0 Qr0 ) ^ wnere we defined = Furthermore, noticing 



that r 



, we can eliminate r and r in favor of fL,, fib, ^r, and £l r in 



(|36|) . Using (fi~8l) . (fl~9~]) and (T20|) we obtain, after some algebra 

7 (z) + A >y( Z y + 4 I 1 



2 V / nr^y 7 (z) 



\ 



(37) 
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where 



with 



i 5i — {1^,0 — n b0 — £l r0 ( 1 + z 



1{Z) -^H Q Q+(z)E*(z) ' [l + zo) ' (38) 



E ( , . _ H (z) = / [(l - <5A0) (1 - O^o - - O-q) + Q b0 ] f l±z_ Y , V r0 fl+z^ 1 



H V 1 - V 1 + z o J 1 - Qj, \1 + z 

(39) 



where A0 (z) = (z) — <po and 5 = A/ /" is an effective coupling constant. As in the 
quintessence field case, if 5 = 0, (157)) reproduces the equation of state parameter obtained 
in Q. 



The evolution of the tachyon scalar field is given by 

d(p _ a/1 + uj^ (z) 



(40) 



dz H E {z){1 + z) 

From (ITS)) . (1T9)) . (120)) and (TTO)) we can calculate the evolution with redshift of all observ- 
ables in the model. If we wish to calculate the time dependence, we need to integrate the 
Friedmann equation (135)) . which can be written in the form 

dt 1 



dz H E(z)(l + z 
From ( )3T)) . we can compute the potential V (z) as 

V(z) 



E 2 (z)Q <t> (z)J-uj^z) , (41) 

PcO v 

where p c0 = SM^H^, E (z) is given by (139)) . (z) is given by (157)) and (z) is the solution 
of flTSJ). From (jHJ) and (j40)), we can compute V{4>). 

The square root in ( 1571) must be real. Furthermore, in analogous manner to the 
quintessence model, must be more than —1 because (jlOj) . We can verify that is 



real and ujs > -1 if (i) < 1 or (ii) > 1 and f- > 2y/E , T a/^^ - 1- 

However, case (ii) is irrelevant, as it corresponds to large values of j^-. For example, if 
fi^o — 0.7 and c = 0.8, we have > 2.45. Below, we will see that the observational data 
constrain ~ 10 _1 . In order that be real for all future times, as — > 1, it is necessary 
that c > 1. As already mentioned in the quintessence case, this is also the condition for the 
entropy to increase for all future times, so the tachyon model also respects the second law 
of thermodynamics. 
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FIG. 3: Equation of state parameter of the holographic tachyon model, for c = 0.85 and = —0.1 
(red dashed line), = (black solid line) and = +0.1 (green dotted line). 

The evolution of the equation of state parameter is showed in figure [3j We have 
ujffi — > —1/3 for z » 1. In the non interacting case -5 = 0- this occurs simply because 
Q(/){z) <C 1 in high redshifts. The behaviour for the interacting case - 5 ^ - is explained as 
follows. In the matter era E 2 (z) ~ (1 + z) 3 so that 7 (z) ~ 1 — -r^-. Using f lTgj) we infer 
^-j^- > 0, that is I7 (2)1 increases with redshift z. Therefore, if 5 < then 7 (z) < and 
increases slower than in the non interacting case. If 5 > then 7 (z) > and increases 
faster than in the non interacting case. In the radiation era 7 (z) ~ an d it turns out to 
be negligible, so that uj^ ~ —1/3. 

In figure HI V(0) is shown for some values of 5 and c. Notice that - as in the case of 
the quintessence field potential - there is a region where V{4>) is almost constant, that is, 
there is a slow-roll region. As we chose <p positive, then evolves to this slow-roll region. 
However, if we had chosen <p negative, then because the right hand side of fj40|) would have 
the opposite sign, so dV ^ would have also the opposite sign and again would evolve to 
the slow-roll region. 

The equation for evolution of <fi ( HU1) can be written in an integral form as 

As before, the model depends on A<p - through E (z) - and neither on <fr nor on <p , then 
it is independent of <po- Therefore, the parameters of the model are 5, c, h, Qbo and Q^o. 
Below, we discuss the comparison with observational data and the results obtained. 
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FIG. 4: Potential of the holographic tachyon field V((f>), in units of p C Q = SM^Hq. cf> is in units 
of Hq 1 . The solid lines are for c = 0.85, the dashed ones are for c = 1.1 and the dotted are for 
c = 1.35. For each value of c the curves from right to left are for = —0.1 (red), = (black) 
and = +0.1 (green), respectively. 

III. CONSTRAINTS FROM OBSERVATIONAL DATA 



In 19j, the lookback time method has been discussed. Given an object i at redshift z i7 
its age t(zi) is denned as the difference between the age of the universe at Z{ and the age of 
the universe at the formation redshift of the object, z F , that is, 

dz' r°° dz' 



t(zi) = H 



-i 



l + z')E(z') 



Ho 1 



(l + z')E(z' 
where ti is the lookback time, given by 



, (l+z')E(z') , Zh 

t L (z F ) - t L (Zi) , 



rr 1 F dz' 

Using (jl2il . the observational lookback time t°^ s {zi) is 



ttiZi) = t L {z F ) ~ t{Zi) = K 8 - t( Zi )} - [t° bS - t L { ZF )\ 



obs 



iobs 



= t° bs -t(zi)-df , 

where is the estimated age of the universe today and df is the delay factor, 



(42) 



(43) 



df = C-t L {z F ) . 
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We now minimize Xmi 



Xm 



N 

£ 

i=i 



[t L (z h p)-tf s ( 



*)] 5 



l 



where ti(zi,p) is the theoretical value of the lookback time in Zi, p denotes the theoretical 
parameters, t°^ s {zi) is the corresponding observational value given by fj43l . <7j is the uncer- 
tainty in the estimated age t(zi) of the object at Zi, which appears in (H3|) and a t obs is the 
uncertainty in getting t^ s . The delay factor df appears because of our ignorance about 
the redshift formation zp of the object and has to be adjusted. Note, however, that the 
theoretical lookback time does not depend on this parameter, and we can marginalize over 
it. 

In [2^] and 21] the ages of 35 and 32 red galaxies are respectively given. For the age 



22( | . Although this estimate for t 



obs 



of the universe one can adopt t^"* = 13.75 ± O.llGyr 
has been obtained assuming a ACDM universe, it does not introduce systematical errors in 
the calculation: any systematical error eventually introduced here would be compensated 
by the adjust of df, in (j43p . On the other hand, such an estimate is in perfect agreement 
with other estimates, which are independent of the cosmological model, as for example 



t° bs = 12.6tl'^Gyr, obtained from globular cluster ages |23j and t^ bs = 12.5 ± 3.0Gyr, 



obtained from radioisotopes studies 24]. 

The WMAP distance information used by the WMAP colaboration includes the "shift 
parameter" R, the "acoustic scale" Ia and the redshift of decoupling z*. These quantities 
are very weakly model dependent [25J]. R and I a are given by 



R = y/Q mO H r(z^) 



and 



I a = rr 



r (2^) 
r s (z*) 



where r (2:*) is the comoving distance to z* and r s (z*) is the comoving sound horizon at z*. 
For a flat universe, r (z*) and r s (z*) are given by 

1 r z * dz 



and 



r s (z*) 



Ho Jo 



Ho Jo E(z) 



dz 



E(z)J3(l + R b /{l + z)) 
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where R b = 3Q b o/ (4Q 7 o). For the redshift of decoupling z* we use the fitting function 
proposed by Hu and Sugiyama [261 ] : 

1 + 0.00124 (Q m h 2 y°™ 



z* = 1048 



+ (n m0 h 2 ) 92 } 



where 



and 



9i 



92 



0.0783 (n b0 h 2 y - 238 
i + 39.5 (n b0 h 2 ) - 763 

0.560 



1 + 21.1 {Q b0 h 2 ) 



1.81 



Thus we add to x 2 the term 



x th — x data ^ 



where x — {I a, R, z*) is the parameter vector and (C *) •• is the inverse covariance matrix 



27J. 



for the seven-year WMAP distance information 

Baryonic Acoustic Oscilations (BAO) are described in terms of the parameter 

2/3 

A= ^iEizBAoY 11 



1/3 


1 


f ZBAO dz' 




_ZBAO 


Jo E{z')\ 



where zbao = 0.35. It has been estimated that A obs = 0.493 ± 0.017 [28]. We thus add to 
X 2 the term 



2 _ {A- ~ A b s ) 
Xbao — ' 



a 



A 



The BAO distance ratio r BAO = D v (z = 0.35) /D v (z = 0.20) = 1.812 + 0.060, estimated 
from the joint analysis of the 2dFGRS (Two Degree Field Galaxy Redshift Survey) and SDSS 
(Sloan Digital Sky Survey) data 29], has also been included. It was demonstrated in 29] 
that this quantity is weakly model dependent. The quantity D v {zbao) is given by 



D v (zbao) = c 
So we have the contribution 



z BAO 



zbao 



\ 2 

A r BAO 



H (zbao) 

_ {tbao - r°i%o) 



dz' 
H(z' 



1/3 



a 



rBAO 



Finally, we add the 397 supernovae data from Constitution compilation 
the distance modulus 

dz' 



301 ] . Defining 



n(z) = 5log 



10 



c(l + z) 



E(z')\ 



+ 25 — 5logioH 
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we have the contribution 



397 



Xsn — 2^ 

3=1 



(J- 



Using the expression x 2 = Xm + Xcmb + Xbao + xl BAO + Xsn, the likelihood function is 



given by 



C(5, c, h^bo^fo) oc exp[ 



X 2 {S, c, h,Q b0 ,tt 



4>o, 



A. Quintessence field 

In table 1 we present the values of the individual best fit parameters, with respective lex, 
2cr and 3<r confidence intervals. 



Table 1: Values of the holographic quintessence model parameters from lookback time, 
CMB, BAO and SNe la. In the last line, Xmm/dof is the minimum \ 2 per degree of 

freedom. 



5M Pl 


n 1 7n+ - 067 + - 148 -+ - 187 

u - 1 ' u -0. 072-0.123-0. 225 


c 


8Q1 +0.048+0.134+0.234 
U.oyi _o.016-0.028-0.035 




fl 77QQ+0.0092+0.0162+0.0264 
V.I '<J'J-0.0087-0.0217-0.0353 




n n/i c;n+ 0026 + 0061 + - 0089 

U.U / ±dU_ .0028-0 .0048-0 .0066 


h 


0.687 ±0.013 ±0.026 ±0.039 


xLJdof 


1.185 



Figure |5] shows the marginalized probability distributions for 5 and c. The coupling 
constant 5 is non vanishing at 2cr confidence level. Figure [6] shows some joint confidence 
regions of two parameters. Notice the effect of the condition < 1 on the positive tail of the 
probability distribuction of 5 and also in the confidence regions of S with other parameters. 
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FIG. 5: Probability distribuctions of the coupling constant 5 (left panel) and of the parameter c 
(right panel) of the holographic quintessence model. 




8M P, c c 

FIG. 6: Two parameters confidence regions of la, 2a and 3a of the holographic quintessence model. 
B. Tachyon scalar field 

In table 2 we present the values of the individual best fit parameters, with respective la, 
2a and 3cr confidence intervals. 
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-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 ' 8 ' g -, ' Q ^ -, ' 2 

5/K r 

o C 

FIG. 7: Probability distribuctions of the coupling constant 5 (left panel) and of the parameter c 
(right panel) of the holographic tachyon model. 

Table 2: Values of the holographic tachyon model parameters from lookback time, CMB, 
BAO and SNe la. In the last line, xLm/dof is the minimum y 2 per degree of freedom. 



s 

Ho 


n 9D1 +0.063+0.117+0.176 
U.ZU1_ .059-0 +38-0 .249 


C 


fl ocq+0.059+0.165+0.235 
U.OOO_ .020-0.030-0.038 




fl 7O++0.013+0.025+0.036 
U. rf ^"±_o. 012-0.026-0. 040 




0.0480 ± 0.0021 ± 0.0041 ± 0.0062 


h 


0.669 ±0.012 ±0.025 ±0.037 


xLJdof 


1.148 



Figure [7] shows the marginalized probability distributions for 5 and c. The coupling 
constant 5 is non vanishing at more than 3a confidence level. So for this model we obtained 
strong evidence for interaction. Figure M shows some confidence regions of two parameters 
for this model. 
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0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 



FIG. 8: Two parameters confidence regions of la, 2a and 3a of the holographic tachyon model. 



IV. DISCUSSION AND CONCLUSIONS 



The minimum \ 2 P er degree of freedom values indicate that both models fit well the 
observational data. Moreover, we see that the tachyon model is a bit more favored by 
the observational data than the quintessence one. The dimensionless coupling constants, 
5Mpi for the quintessence field and -J- for the tachyon, agree at la level and for both 
models we obtained significative evidence for a non vanishing interaction in the dark sector. 
Furthermore, both the results implies in dark energy decaying into dark matter, alleviating 
the coincidence problem. These results are consistent with previous ones, as for example 
those obtained in {4], where an interacting tachyonic dark energy model with a power law 
potential was assumed and the results were consistent with the interaction at 90% confidence 
level. We must also mention that in Q and evidence for interaction was found using 
completely different models and data sets. So those results combined with the present ones 
furnish sensible evidence in favor of an interaction in the dark sector of the universe. 

The results obtained for c, O^o and h for the two models also agree at la confidence 
level. The values obtained for f^ agree only at 2a level. For the quintessence model, fl^o is 
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almost superestimated, corresponding to a matter relative density today O m o = 0.231±0.010. 
This value agrees at inferior limit of la with the cosmological model independent estimative 
^Mobs = 0.28±0.04 [3l|. For the tachyon model, Q^q corresponds to a matter relative density 
today Q m o = 0.276^g'o^25 which is in perfect agreement with that observational estimative. 
For both models, the baryionic density and the Hubble parameter today are very reasonable. 
We have obtained Q b0 h 2 = 0.0212±0.0015 and Q b0 h 2 = 0.0215±0.0012 from quintessence and 
tachyon models, respectively. Let's compare these values, for example, with that obtained 
from deuterium to hydrogen abundance ratio [32|, flboh 2 = 0.0213 ± 0.0013 ± 0.0004, where 
the errors terms represents the la errors from deuterium to hydrogen abundance ratio and 
the uncertainties in the nuclear reaction rates, respectively. For the Hubble parameter, we 
have obtained h = 0.687 ± 0.013 and h = 0.669 ± 0.012 from the quintessence and tachyon 
models, respectively, both in excellent agreement with observational values, independent 



of cosmological model, as for example h Q b s = 0.69 ± 0.12 [20j and h b s = 0.72 ± 0.08 [33 ]. 
We can also compare the ratios fi&o/^mo = 0.248^q'q23, from the quintessence model, and 
^f>o/^mo = 0.211 ± 0.016, from the tachyon model, with the observational value of 2dFGRS 
colaboration, Q bo /Q m o = 0.185 ± 0.046 [34|. 

We have obtained c < 1 at la confidence level for both models. As already said above, 
this implies that the equation of state parameter will not be real for all future times. 
However, this is not a very serious problem, because c is compatible with values above unit at 
2a confidence level. Moreover, one could say that c < 1 is only an effect due to lack of more 
precise observational data. Anyway, the very simple models presented here are expected 
to be only alternatives to an effective description of a more sophisticated subjacent theory 
of dark energy. In principle, nothing guarantees that they will be good descriptions for all 
future times. 

Figures [6] and [8] shows some joint confidence regions of two parameters for both models. 
In the confidence regions for 5 versus c and for c versus Q<f>o, we see that there is a lower 
limit on c ~ 0.8. This also can be seen in the marginalized probability distributions of c, 
which dies for c < 0.85. This lower limit is explained by the condition < 1, necessary 

for ojfi to be real and > —1, discussed above. This limit can be seen more clearly in c 
versus f2^ confidence regions. Moreover, we have c ~ a/^o f° r the best fit values of these 
parameters. This implies that u^o — —1 and both models approaches AC DM today. This 
is consistent with the fact that, as ACDM fits all observational data, then any alternative 
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model must not deviates much from ACDM for z ~ 0. However, for z > 0, both models 
are qualitatively different from ACDM. For the quintessence field we have very different 
qualitative behaviours for 5 < 0, 5 = and 5 > 0, as showed in figure [T] and discussed in 
11 A. For the tachyon model its behaviour is qualitatively the same in the three cases, 
approximates —1/3, as showed in figure [3] and discussed in II B. 

It is interesting to compare the results obtained in the present work for the holographic 
tachyon model with the previous ones, presented in [6|, where a simpler version of the 
holographic tachyon model, without barions nor radiation, had been compared with obser- 
vational data. The values obtained here for c, VL^q and h are the same as before and with 
minor incertainty intervals, despite the fact that here we have one more parameter - -> 
as we can see by comparing table 2 in the present work with table 1 in jf}]. This is because 
now it was possible to use all WMAP distance information R, I a and z*, as the model was 
generalized to include barions and radiation. However the dark energy coupling constant 5 
now is non vanishing and compatible with dark energy decaying into dark matter with more 
than 3cr confidence level, whereas in 6|] no evidence of interaction had been found. This 
can be understood as follows. In low redshifts, the universe is dominated by dark energy. 
The dynamics of dark energy in low redshifts is essentially determined by the equation of 
state parameter ^(z), as can be seen in (fT8l) . As in this period the data sets were the 
same in both the works, so (z) must be almost the same in this period in both the works. 
But Co>0 (z) explicitly depends on the product SQ^o, were 6 is the coupling constant and 
Qq,o = 1 — Q^o — &>bo — ^ro is the dark matter relative density, as we can see in ( 1371) and ( 1381) . 
As in the present work Q^q is less than that in [6] - where flu) = Q r o = -, it turns out that 
5 in the present work is bigger (in modulus) than that in 6j. However it is important to 
point out that in [g| the model was less realistic, as it wasn't include barions nor radiation. 
Therefore, the present result favorable to interaction is more robust. 

In summary, combinations of holographic dark energy model and scalar fields were im- 
plemented. It was showed that it is possible to fix the potential of interacting scalar fields 
by imposing that the energy density of the scalar field must match the energy density of 
the holographic dark energy. A comparison of the models with recent observational data 
was made and the coupling is non vanishing at more than 2a for the quintessence field and 
at more than 3a for the tachyon. In both cases the results are consistent with dark energy 
decaying into dark matter, alleviating the coincidence problem. 
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